Enumeration of many-body skeleton diagrams 
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The many-body dynamics of interacting electrons in condensed matter and quantum chemistry 
is often studied at the quasiparticle level, where the perturbative diagrammatic series is partially 
resummed. Based on Hedin's equations for self-energy, polarization, propagator, effective potential, 
and vertex function in zero dimension of space-time, dressed Feynman (skeleton) diagrams are 
enumerated. Such diagram counts provide useful basic checks for extensions of the theory for future 
realistic simulations. 
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I. INTRODUCTION 

Current research in electronic properties of molecules, 
nanocomponents, and solids has gone far beyond the 
mean-field density-functional description. Many-body 
methods are employed routinely at the level of quasi- 
particles, to describe excitations accuratelyi^. Propa- 
gators and interactions are renormalized to obtain an ef- 
fective theory of weakly interacting quasiparticles. This 
is achieved by infinite resummations of Feynman dia- 
grams, that are brought to a smaller class of skeleton 
diagrams 4 . The reduction in number should correspond 
to better analytic behavior in space-time of the individ- 
ual dressed diagrams, as it occurs in random-phase or 
ladder resummations^. 

In this report, we show how skeleton Feynman dia- 
grams of various types can be enumerated for the exact 
many-body theory and its GW approximation 1 * 2 ^. We 
take as our starting point the standard formulation of the 
many-body problem of TV interacting fermions provided 
by the set of five Hedin's equationsSiiSiX for the propa- 
gator G(l, 2), effective potential W(l, 2), irreducible self- 
energy £(1,2), irreducible polarization 11(1,2), and irre- 
ducible vertex T(l; 2, 3). 

In zero dimension of space and time, the combina- 
torial content of Wick's expansion survives, producing 
the same Feynman diagrams in the perturbative expan- 
sion of the correlators as in conventional 3+1 dimensions. 
These correlators no longer carry space-time labels, and 




FIG. 1: The unique fully dressed skeleton diagrams for self- 
energy (left) and polarization (right). 



do not correspond to space-time functions. However, 
they continue to solve Hedin's equations, which can in- 
deed be derived from considerations about the topology 
of diagrams^. In zero dimensions (d = 0), four of Hedin's 
equations become algebraic, in terms of the scalar vari- 
ables g and v, representing the Hartree propagator (com- 
puted with the exact density provided by G) and the 
bare interparticle interaction respectively. The functional 
derivative in the equation for the vertex becomes an or- 
dinary derivative. The five equations read: 

G = g + gY,G, W = v + vUW , (1) 

s =iGwv, n = ^G 2 r, (2) 
r =i + rG 2 §§. (3) 

The parameter £ is introduced to count fermion loops: 
in Hedin's equations for the electron, £ = — 2 because 
of spin degeneracy and Feynman's rule for the fermionic 
loop. The diagrammatic meaning of Hedin's equations is 
simple: Eqs. Q correspond to Dyson's equations that de- 
fine the proper self-energy and polarization insertions for 
the propagator and the effective potential, Eqs. (0) trans- 
late the unique skeleton structure^ of £ and n, shown in 
Fig. Finally, the vertex equation © shows that vertex 
diagrams arise from self-energy diagrams, by "pinching" 
a G line with a vertex (remove a line G and replace it 
with GrG): this follows from the functional definition 
of the vertei&&&. In the GW approximation all correc- 
tions to the bare vertex are neglected: the approximation 
r = 1 replaces the exact vertex equation (0. 

II. ENUMERATION OF FEYNMAN 
DIAGRAMS 

The perturbative solution of Hedin's equations in zero 
dimension provides numerical coefficients which enumer- 
ate the Feynman diagrams for the five correlators in- 
volved. The ordinary perturbation expansion is carried 
out in the bare parameters v and g. However, it is often 
convenient to expand in different "renormalized" vari- 
ables, such as v and G, or W and G, or G,W and T. These 
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TABLE I: The number of °tl -order skeleton diagrams for the vertex function in the five renormalization schemes considered 
in the text. 
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r(x) 
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342430 
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63 


729 


10113 


161935 


2923135 
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30747331223 


r(t) 
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6 


52 


602 


8223 


128917 


2273716 


44509914 
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22449011336 


r(«) 


i 


6 


49 


542 


7278 


113824 


2017881 


39842934 


865391422 


20486717908 


i» 
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3 


13 


147 


1965 


30979 


559357 


11289219 


250794109 


6066778627 




FIG. 2: The four diagrams for G at second order in x. 



The GW approximation (T = 1), generates a cubic equa- 
tion for Zgw- For the complete theory, we rewrite the 
vertex equation Q, in a form^ where G is traded for g: 

^ ^ ^ dg Z{x) dx Z(x) ^ 

Equations (@J and JHJ combine to a closed equation for 
Z: 



expansions count Feynman diagrams where, respectively, 
the propagator or both the propagator and the potential, 
or also the vertex, have been fully rcnormalizcd. Such di- 
agrams where self-energy contributions of various types 
are resummed, are called skeleton graphs^. Different 
levels of resummations will be considered below, includ- 
ing a renormalization based on GW correlators. 
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with the initial condition Z(0) = 1. The solution as 
a series expansion in x provides the number of Feynman 
diagrams that contribute to each order to the propagator 
G 



x = g u-expansion 



Z(x) = 1 + x + (3 + £)x 2 + (15 + lU + £ 2 )x 3 



(7) 



Standard perturbation theory in the bare interaction 
v has been studied in d = by various authorsi&Ai by 
considering the path integral formulation of the interact- 
ing field, which degenerates to an ordinary integral. An 
approach based on Hedin's equations is very convenient 
to deal with the Hartree propagator g in place of the bare 
one, and explicit counting numbers can be obtained^. 

Here, the natural dimensionless expansion parameter 
is x = ig 2 v. As an example, we provide the equation for 
Z(x) — G/g, where Z is the renormalization factor due 
to self-energy corrections. The algebraic equations (11121) 
yield 



lxZ 2 {x)T{x) 



1 



1 



1-1 + iZ(x) 



(4) 



The £ powers represent the loop contents. For example, 
at second order in x, the theory involves 3 diagrams with 
no loops, and 1 diagram with one loop; they are shown 
in Fig. [3 By substituting the power series for Z(x) into 
Eq. © , we obtain the diagram count for the vertex func- 
tion: 



T{x) = l + x + 3(2 + £) x 2 + 5(l0 + 9£ + £ 2 )x 3 + ...(8) 

Table [I] reports the total number of vertex diagrams up 
to 10 th order, obtained by taking £ = 1 in the expansion 
for T(x). 



B. y — G 2 u-skeleton expansion 

When all self-energy insertions of the propagator are 
resummed, one ends up with Feynman diagrams where 



the propagator lines correspond to the exact G, while 
the interaction lines continue to correspond to the bare 
v. This resummation is obtained by using the expansion 
parameter y — iG 2 v. It is possible to obtain a closed 
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equation for the vertex function. First write the self- 
energy as 



£ = iGWT = %G 



-r = 



i y r 



l-vll Gl-ylT 



(9) 



Next evaluate the vertex function in Eq. J3J by means of 
Eq. 10, and obtain a differential equation for T(y): 



2y 2 T 



dy 



-1 + (1 + 2y£)T - yr 2 (l + 2£ + y£ 2 ) 
+y 2 T 3 e(£-l). (10) 



This equation must be solved with the initial condition 
r(0) = 1. The series expansion of the solution counts 
skeleton vertex diagrams with dressed propagators and 
bare interactions: 

T(y) = 1 + y + (4 + 3£)y 2 + (27 + 3U + U 2 )y 3 + . . (11) 

Table [I] lists the total number of these diagrams up to 
order n = 10, by taking I = 1 in the expansion for r(y), 
as was done above. 

The enumeration of skeleton diagrams of the polar- 
ization II = i£G 2 T(y) coincides with that of T(y). The 
skeleton expansion of £ results from Eq. © : 

E/iGv = 1 + (1 + £)y + (4 + U + £ 2 )y 2 

+ (27 + 40^+14£ 2 +^ 3 )y 3 + ... (12) 

In the expansion in terms of the bare parameter x, H 
and £ shared the diagram counting,- while this property 
does not apply to the y expansion at hand. 

The G 2 v-expansion of £ is interesting for the theory 
of the Luttinger-Ward ^-functionalist. By closing all 
£[G, v] skeleton graphs with a G line, and dividing each 
one by the number of G lines it contains, one obtains a 
functional with the property of yielding the self-energy 
and the reducible polarization: 

In d = 0, the functional ^[G, v] turns into an ordinary 
function of y = iG 2 v. 
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2j ~ a 1 '//T: /; 
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- f + (l+^ + (4 + 4^ 2 )^ + ... (14) 

This last expression provides the appropriate diagram 
counting and fractional weights for $(?/) in any dimen- 
sion. 



C. z — G 2 W-skeleton expansion 

Resummation of both self-energy and polarization in- 
sertions leads to skeleton diagrams with exact propaga- 
tors G and interactions W. Here, the natural expansion 



parameter is z = iG 2 W. The parameter z is linked as 
follows 



iG 2 



y 



1 - i£vTG 2 1 - £yT 



(15) 



to the parameter y — iG 2 v of the previous expansion 
(|llfl . Inversion yields y — z [1 + £zT(z)]~ 1 . By entering 
this relation into the differential equation (|10(l for T(y), 
which is now viewed as a function of z, we obtain 

z2 dT(z) = l-T(z) + zT 2 (z) + 2£z 2 r 3 (z) 



dz 



£-{2 + £)T(z)~iz£V 2 {z) 



The appropriate series expansion counts G 2 M^-dressed 
vertex diagrams: 



T(z) = 1 + z + 2 (2 + £) z 2 + (27 + 22£) z 3 + 



(17) 



Table [I] lists the total number of these T(z) diagrams up 
to order n = 10. Both £ = iGWT(z) and II = ^G 2 r(z) 
have the same total G 2 Vy-counting numbers as the ver- 
tex. 

In the G 2 II 7 -expansion, it is natural to define^ the 
functional ^[G, W), which generalizes Eq. ((T^l and gen- 
erates the self-energy and the irreducible polarization as 
follows: 

S(1 ' 2) ~ JG(2,1) ' n(1 ' 2) -" 2 W(2,l) • (18) 

The diagrammatic construction of ^ is the same as for $: 
close the skeleton graphs of £ = iGWT[G, W] with a G 
line and divide by the number of G lines that the graph 
contains. In d — we have G£ = zT(z); to divide by 
the number of G lines that each G£-diagram contains, 
corresponds to the integral 



Viz) 



1 



D. 



2 3 

dz'T(z') ^^+ 7 — + {2 + £) Z — + ... (19) 



u = G 2 WT 2 -skeleton expansion 



The G 2 W— expansion resums all self-energy and polar- 
ization insertions, and skeleton diagrams of various order 
in z result just because of vertex contributions. If vertex 
diagrams are summed as well, each of the self-energy and 
the polarization is brought to a unique skeleton diagram 
depicted in Fig. ^ Vertex diagrams cannot be brought 
to a finite collection of skeleton diagrams. We may nev- 
ertheless resum vertex insertions in all vertex diagrams, 
and enumerate the resulting vertex skeleton diagrams, 
where all lines and vertices are resummed. 

This is achieved by noting that in a vertex diagram, 
each W line ends in two vertices. However, the vertex 
which the Coulomb external line is attached to, is left 
out (see Fig. [3] for some vertex diagrams). We then write 
r = 1 + T"/(u), in terms of the expansion variable u = 
iG 2 WT 2 — zT 2 (z). The factor V in front of j(u) accounts 
for the left-out vertex. 
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FIG. 3: All first-order (a) and second-order (b,c,d) vertex 
diagrams in u expansion. 



The Taylor expansion of j(u) enumerates all dressed 
fully renormalized vertex diagrams. An equation for 7 is 
obtained from the relation 



(hi 
du 



dz d r — 1 



1 



dT I dz 



dudz T T 2 T 2 + 2zTdT/dz 
where we enter dT/dz given by Eq. Ijltjfl : 



dj 
du 



(1-7) 



2£u 2 (l-7) 2 + u(l-7) -7 
£u 2 (l - 7 ) 2 - £wy{l - 7) - 2 7 



(20) 



(21) 



This equation is solved with initial condition 7(0) = 0. 
The resulting u-expansion for T is 



r(«) 



1 + (Ttt) + (1 + 2£)(Tu 2 ) + (7 + 6£)(Tu 3 ) 
+(63 + 7A£ + 10£ 2 ){Tu 4 ) + ... (22) 



The diagrams of first and second order are shown in Fig. [3] 
(the next 13 third-order diagrams coincide with those of 
QED, drawn in Ref. llOt) . This u-expansion represents 
the "ultimate" skeleton expansion, where all ingredients 
of Hedin's equations have been renormalized. Indeed, as 
apparent in Table [IJ the diagram count is smallest in the 
u expansion at hand. 



E. t = (G 2 W / ) G w-skeleton expansion 

In the GW approximation (r = 1), in physical dimen- 
sions, Hedin's equations are a system of ordinary integral 
equations for Gqw, ^gw, Wgw and IIgh/. It is conceiv- 
able that in a near future increased computer power will 
allow us to solve these equations in fully self-consistent 
GW: this would make the GW approximation the zeroth- 
order stage of a subsequent attack of the full many-body 
problem, where vertex corrections are included pertur- 
batively. The problem of including vertex corrections 



03 0.6 



0.4 



3 
9 
a 



0.2 



"1 1 I ' T 
x-expansion 



y-expansion 



t-expansion 
»« >>♦ « < 




z-expansion 
u-expansion 



_L 



_L 



10 



20 30 
order n 



40 



50 



FIG. 4: (color online) Number of vertex diagrams at order n, 
divided by the asymptotic expression B n defined in Eq. 1271 , 
for the different renormalization schemes realized in Sec. Ill Al 
[x - Eq. 0], Sec.|TTU[y - Eq. JJTJ], Sec.|TTEl[f - Eq. JUJ], 
Sec. Hill [z - Eq. C0], and Sec. HTT31 [u - Eq. JUJ]. 



in a systematic way is an important target of present- 
day research, to overcome several limitations of the non- 
conservative GWii^i^. It is therefore of interest to count 
the diagrams where g and v lines are dressed (resummed) 
to include all their GW self-energy and polarization in- 
sertions. 

The counting problem is solved in d = by considering 
the expansion parameter t = iG 2 GW Wcw- With simple 
algebra, the following two equations for the vertex and 
the ratio Zcw(t) = G/Ggw are obtained: 



tz, 



,;w(t)m 
dZcw 



2f~ 



dt 



[(l + t)Z GW {t)-l](l+£t) 
I -£ + £{! + t)Z Gw it) 

{ l-t-2£t 2 ) 



(23) 



t 



l + it 



\-£ + £{l + t)Z GW (t) 
Z G w(t). (24) 



The perturbative solutions are: 

Z G w(t) = 1 + t 2 + {5 + 3£)t 3 + ... (25) 
T(t) = l + t + 2{2 + £)t 2 + (29 + 23£)t 3 + . .(26) 

Table lists the total number of these T(t) diagrams up 
to order n — 10. 



III. ASYMPTOTICS 

To estimate asymptotic behaviors is not simple in the 
present approach. However, the topology of many-body 
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diagrams with two-body interaction is the same as in 
relativistic quantum electrodynamics (QED), with a dif- 
ference: because of exact particle-antiparticle symme- 
try, QED diagrams with loops involving an odd number 
of fermion lines cancel (Furry's theorem^). Therefore, 
many-body skeleton diagrams grow faster in number than 
QED ones. 

The functional approach, in the saddle point expan- 
sion, is better suited for asymptotics. This problem in 
physical dimensions has been studied for QED by sev- 
eral authors, based on Lipatov's powerful methodiSii!. 
From the tabulated QED asymptotic values^ in d = 0, 
one can infer the leading behavior of many-body ver- 
tex skeleton diagrams. Consider the series for the vertex 
= Xin^™^™) m one °f the renormalization schemes 
(9 = x, y, z, t, u) outlined above. We obtain 

A n = B n (c + c x /n +...), B n = n\ 2" n 3 / 2 . (27) 

In Fig. 0] the asymptotic behavior is checked on the first 
50 coefficients of the five different vertex skeleton expan- 
sions computed in the present work. Indeed, in all ex- 
pansions, whether bare or renormalized, the number of 
vertex diagrams grows with the order n roughly as B n , 
thus leading to divergent power series, like in QED. The 
renormalization scheme only affects the subleading co- 
efficients Cj. As expected, the "ultimate" u expansion, 
involving wider diagram resummations leads to fewer dia- 



grams than all other "less renormalized" expansions. For 
large n, the diagram count is essentially the same for 
the z expansion, including renormalization to both the 
propagator and the interaction, and for the t expansion 
around Hedin's GW approximation. 



IV. CONCLUSIONS 

The present method for enumerating skeleton diagrams 
is based on Hedin's equations and it is very efficient. 
Simple changes of the expansion variable produce var- 
ious renormalizations schemes, in the form of differential 
equations that can easily be solved by series. The coeffi- 
cients of these series are integers that count skeleton dia- 
grams. These are useful as a check for any diagrammatic 
approach to the many-body problem in realistic dimen- 
sions. As expected, the number of diagrams, whether 
renormalized or not, grows factorially with the order. 
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